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Abstract Strongly interacting Fermi gases provide a clean and controllable lab- 
oratory system for modeling strong interparticle interactions between fermions in 
nature, from high temperature superconductors to neutron matter and quark-gluon 
plasmas. Model-independent thermodynamic measurements, which do not require 
theoretical models for calibrations, are very important for exploring this impor- 
tant system experimentally, as they enable direct tests of predictions based on the 
best current non-perturbative many-body theories. At Duke University, we use 
all-optical methods to produce a strongly interacting Fermi gas of spin-l/2-up and 
spin-l/2-down 6 Li atoms that is magnetically tuned near a collisional (Feshbach) 
resonance. We conduct a series of measurements on the thermodynamic properties 
of this unique quantum gas, including the energy E, entropy S, and sound veloc- 
ity c. Our model-independent measurements of E and S enable a precision study 
of the finite temperature thermodynamics. The E(S) data are directly compared to 
several recent predictions. The temperature in both the superfluid and normal fluid 
regime is obtained from the fundamental thermodynamic relation T = dE/dS by 
parameterizing the E(S) data using two different power laws that are joined with 
continuous E and T at a certain entropy S c , where the fit is optimized. We ob- 
serve a significant change in the scaling of E with S above and below S c . Taking 
the fitted value of S c as an estimate of the critical entropy for a superfluid-normal 
fluid phase transition in the strongly interacting Fermi gas, we estimate the critical 
parameters. Our E(S) data are also used to experimentally calibrate the endpoint 
temperatures obtained for adiabatic sweeps of the magnetic field between the ideal 
and strongly interacting regimes. This enables the first experimental calibration of 
the temperature scale used in experiments on fermionic pair condensation, where 
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the ideal Fermi gas temperature is measured before sweeping the magnetic field 
to the strongly interacting regime. Our calibration shows that the ideal gas tem- 
perature measured for the onset of pair condensation corresponds closely to the 
critical temperature T c estimated in the strongly interacting regime from the fits 
to our E(S) data. We also calibrate the empirical temperature employed in stud- 
ies of the heat capacity and obtain nearly the same T c . We determine the ground 
state energy by three different methods, using sound velocity measurements, by 
extrapolating E(S) to S = and by measuring the ratio of the cloud sizes in the 
strongly and weakly interacting regimes. The results are in very good agreement 
with recent predictions. Finally, using universal thermodynamic relations, we es- 
timate the chemical potential and heat capacity of the trapped gas from the E(S) 
data. 

PACS numbers: 03.75.Ss 



1 Introduction 

Interacting fermionic particles play a central role in the structure of matter and 
exist over a very broad range of energies, from extremely low temperature trapped 
atomic Fermi gases, where T < 10~ 7 K^* 2 -, to very high temperature primordial 
matter, like quark-gluon plasmas, where T > 10 12 K— , For all of these systems, 
the most intriguing physics is related to very strong interactions between fermionic 
particles, such as the strong coupling between electrons in high-T c superconduc- 
tors and the strong interactions between neutrons in neutron matter. 

Current many-body quantum theories face great challenges in solving prob- 
lems for strongly interacting Fermi systems, due to the lack of a small coupling 
parameter. For example, the critical temperature of a superfluid-normal fluid tran- 
sition in a strongly interacting Fermi gas has been controversial for many years. 
The critical temperature T c /Tp has been predicted to have values in the range be- 
tween 0.15 and 0.35 by different theoretical method s 4 i 5 1 6 i 7 1 8 i 9 1 1 ° 1 1 1 . A complete 
understanding of the physics of strongly interacting systems can not yet be ob- 
tained from a theoretical point of view. There is a pressing need to investigate 
strongly interacting fermions experimentally. 

In recent years, based on progress in optical cooling and trapping of fermionic 
atoms, a clean and controllable strongly interacting Fermi system, comprising 
a degenerate, strongly interacting Fermi gas 1 ' 2 , is now of interest to the whole 
physics community. Strongly interacting Fermi gases are produced near a broad 
Feshbach resonance W 3 , where the zero energy s-wave scattering length as is 
large compared to the interparticle spacing, while the interparticle spacing is large 
compared to the range of the two-body interaction. In this regime, the system is 
known as a unitary Fermi gas, where the properties are universal and indepen- 
dent of the details of the two-body scattering interactio n 14 i 15 . In contrast to other 
strongly interacting Fermi systems, in atomic gases, the interactions, energy, and 
spin population can be precisely adjusted, enabling a variety of experiments for 
exploring this model system. 

Intense studies of strongly interacting Fermi gases have been implemented 
over the past several years from a variety of perspectives. Some of the first experi- 
ments observed the expansion hydrodynamics of the strongly interacting cloud^*^. 
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Evidence for superfluid hydrodynamics was first observed in collective mode s ' . 
Collective modes were later used to study the T = equation of state throughout 
the crossover regim e 19 i 20 i 21 . Recently, measurements of sound velocity have also 
been used to explore the T = equation of state 22 . Below a Feshbach resonance, 
fermionic atoms join to form stable molecules and molecular Bose-Einstein con- 
densate s 23 i 24 i 25 1 26 i 27 . Fermionic pair condensation has been observed by projec- 
tion experiments using fast magnetic field sweep s 26 i 27 . Above resonance, strongly 
bound pairs have been probed by radio frequency and optical spectroscop y 28 i 29 i 30 i 3 
Phase separation has been observed in spin polarized samples 32 ' 33 . Rotating Fermi 
gases have revealed vortex lattices in the superfluid regim e 34 ' 35 as well as irrota- 
tional flow in both the superfluid and normal fluid regimes 36 . Measurement of the 
thermodynamic properties of a strongly interacting Fermi gas was first accom- 
plished by adding a known energy to the gas, and then determining an empiri- 
cal temperature that was calibrated using a pseudogap theory 37 . Recent model- 
independent measurements of the energy and entropy 38 provide a very important 
piece of the puzzle, because they enable direct and precision tests that distinguish 
predictions from recent many-body theories, without invoking any specific theo- 
retical model^. 

One of the major challenges for the experiments in strongly interacting Fermi 
gases is the lack of a precise model-independent thermometry. Two widely-used 
thermometry methods are model-dependent, in that they rely on theoretical mod- 
els for calibration. The first relies on adiabatic magnetic field sweeps between the 
molecular BEC regime and the strongly interacting regime 28 - 39 . Subsequently, the 
temperature of the strongly interacting gas is estimated from the measured tem- 
perature in the BEC regime using a theoretical model of the entropy 8 . The second 
method, used by our group 37 , is based on determining an empirical temperature 
from the cloud profiles that is calibrated by comparing the measured density dis- 
tribution with a theoretical model for the density profiles. Currently two model- 
independent thermometry methods have been reported for strongly-interacting 
gases. One is the technique employed by the MIT group 40 , which is only applica- 
ble to imbalanced mixtures of spin-up and spin-down atoms. That method is based 
on fitting the noninteracting edge for the majority spin after phase separation. An- 
other model-independent method is demonstrated in Ref. 38 , which is applicable 
to both balanced and imbalanced mixtures of spin-up and spin-down fermions. 
The energy E and entropy S are measured and then parameterized to determine a 
smooth curve E(S). Then the temperature in both the superfluid and normal fluid 
regime is obtained from the fundamental thermodynamic relation T = dE/dS. 

In this paper, we will describe our model-independent thermodynamic exper- 
iments on a strongly interacting Fermi gas of 6 Li, which we have conducted at 
Duke University. First, we will describe our measurements of both the total en- 
ergy E and the total entropy S of a trapped strongly-interacting Fermi gas tuned 
near a Feshbach resonance. Then, we determine the temperature T = dE/dS af- 
ter showing that the E(S) data are very well parameterized by using two different 
power laws that are joined with continuous E and T at a certain entropy S c that 
gives the best fit. To examine the sensitivity of the temperature to the form of the 
fit function, we employ two different fit functions that allow for a heat capacity 
jump or for a continuous heat capacity at S c . We find that the T values closely 
agree for both cases. We find a significant change in the scaling of E with S above 
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and below S c , in contrast to the behavior for an ideal Fermi gas, where a single 
power-law well parameterizes E(S) over the same energy range. By interpreting 
S c as the critical entropy for a superfiuid-normal fluid transition in the strongly 
interacting Fermi gas, we estimate the critical energy E c and critical temperature 
T c . Both the model-independent E(S) data and the estimated critical parameters 
are compared with several recent many-body theories based on both analytic and 
quantum Monte Carlo methods. 

We also show how parameterizing the E(S) data provides experimental tem- 
perature calibrations, which helps to unify, in a model-independent way, the results 
obtained by several groups 26 ,27 i 31 ' 37 ' 38 . First we relate the endpoint temperatures 
for adiabatic sweeps of the bias magnetic field between the strongly interacting 
and ideal noninteracting regimes, as used in the JILA experiments to characterize 
the condensed pair fractio n 26 ' 31 . This enables the ideal gas temperature observed 
for the onset of pair condensation 26,31 to be related to the critical temperature of 
the strongly interacting Fermi gas. The temperature obtained by parameterizing 
the strongly interacting gas data also calibrates the empirical temperature based 
on the cloud profiles, as used in our previous studies of the heat capacity 37 . These 
temperature calibrations yield values of T c close to that estimated from our E(S) 
data. 

Next, we discuss three different methods for determining the universal many- 
body parameter, j3 1 , where 1 + j3 is the energy per particle in a uniform strongly 
interacting Fermi gas at T = in units of the energy per particle of an ideal Fermi 
gas at the same density. First, we describe the measurement of the sound velocity 
at resonance and its relationship to j3. Then, we determine j3 from the ground 
state energy Eq of the trapped gas. Here, Eo is obtained by extrapolating the E(S) 
data to S = 0, as suggested by Hu et al.— . This avoids a systematic error in the 
sound velocity experiments arising from the unknown finite temperature. Finally, 
to explore the systematic error arising from the measurement of the number of 
atoms, P is determined in a number-independent manner from the ratio of the 
cloud sizes in the strongly and weakly interacting regimes. All three results are 
found to be in very good agreement with each other and with recent predictions. 

Finally, we obtain three universal thermodynamic functions from the param- 
eterized E(S) data, the energy E(T), heat capacity C(T), and global chemical 
potential n g (E). 



2 Experimental Methods 

Our experiments begin with an optically-trapped highly degenerate, strongly in- 
teracting Fermi gas of 6 Li 1 . A 50:50 mixture of the two lowest hyperfine states of 
6 Li atoms is confined in an ultrastable CO2 laser trap with a bias magnetic field 
of 840 G, just above a broad Feshbach resonance at B = 834 G^l. At 840 G, the 
gas is cooled close to the ground state by lowering the trap depth [7 1,13 . Then 
U is recompressed to a final trap depth of Uq/Icb = 10;UK, which is much larger 
than the energy per particle of the gas, for the highest energies employed in the 
experiments. This suppresses evaporation during the time scale of the measure- 
ments. The shallow trap yields a low density that suppresses three body loss and 
heating. The low density also yields a weakly interacting sample when the bias 
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magnetic field is swept to 1200 G, although the scattering length is —2900 bohr, 
as discussed in detail in § 12.21 

The shape of the trapping potential is that of a gaussian laser beam, with a 
transverse gaussian profile determined by the spot size and an axial lorentzian 
profile determined by the Rayleigh length. To simplify the calculations of the ideal 
gas properties in subsequent sections, as well as the theoretical modelling, we take 
the trap potential to be approximated by a three dimensional gaussian profile, 



where a Xl y :Z is the 1/e 2 width of trap for each direction. Here, we take the zero of 
energy to be at r = 0. When the cold atoms stay in the deepest portion of the optical 
trap, where x(y, z) << a x (a y ,a z ), the gaussian potential can be well approximated 
as a harmonic trap with transverse frequencies CO x , COy and axial frequency ft),, 
where 



Here m is the 6 Li mass. At our final trap depth Uo, the measured transverse fre- 
quencies are co x = 2nx 665(2) Hz and co y — 2% x 764(2) Hz. The axial frequency 
is weakly magnetic field dependent since the total axial frequency has both an 
optical potential contribution ft) 2 ., determined by Eq. [2] and a magnetic poten- 
tial contribution arising from magnetic field curvature, ft) 2 ir The net axial fre- 
quency is then ft) z = \J + 0)} nz . We find ft), = 2iz x 30.1(0.1) Hz at 840 G and 
ft) z = 2n x 33.2(0.1) Hz at 1200 G. The total number of atoms is N ^ 1.3 x 10 5 . 
The corresponding Fermi energy Ep and Fermi temperature Tp at the trap center 
for an ideal noninteracting harmonically trapped gas are Zip = kgTp = hcd(3N) 1 ' 3 , 
where ft) = (ft^ft^ft),) 1 / 3 . For our trap conditions, we obtain Tp ~ 1.0/lK. 
Using ft), we can rewrite Eq.Q]as a symmetric effective potential, 



where r is the scaled position vector. Here, r 2 = x 2 +y 2 +z 2 with* = co x x/co,y = 
(Q y y/(D,z= CO z z/cd. To obtain the anharmonic corrections for the gaussian trap, we 
expand Eq.[3]in a Taylor series up to second order in r 2 , 
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2. 1 Energy Measurement 



Model-independent energy measurement is based on a virial theorem, which for 
an ideal gas in a harmonic confining potential Uho yields E = 2(Uho)- Since the 
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harmonic potential energy is proportional to the mean square cloud size, measure- 
ment of the cloud profile determines the total energy. Remarkably, a trapped uni- 
tary Fermi gas at a broad Feshbach resonance obeys the same virial theorem as an 
ideal gas, although it contains superfluid pairs, noncondensed pairs, and unpaired 
atoms, all strongly interacting. This has been demonstrated both theoretically and 
experimentally 42 . The virial theorem shows that the total energy of the gas at all 
temperatures can be measured from the cloud profile using 



where U is the trapping potential and x is the position vector. Eq. [5]can be shown 
to be valid for any trapping potential U and for any spin mixture, without assuming 
either the local density approximation or harmonic confinemen t 43 ' 44 ' 45 ' 46 . 

Using Eq.|4]in Eq.[5]and keeping the lowest order anharmonic corrections, we 
obtain the energy per atom in terms of the axial mean square size, 



Here, we have used the local density approximation with a scalar pressure, which 
ensures that (xdU /dx) = (ydU /dy) = (zdU /dz). For the ground state, where 
the spatial profile is a zero temperature Thomas-Fermi profile, we have (z 4 ) = 
12(z 2 ) 2 /5. For energies EjEp > 1, where the spatial profile is approximately 
gaussian, we have (z 4 ) = 3(z 2 ) 2 . Since the anharmonic correction is small at low 
temperatures where the cloud size is small, we use the gaussian approximation 
over the whole range of energies explored in our experiments. For the condi- 
tions of our experiments, there is no evidence that the local density approxima- 
tion breaks down for a 50:50 spin mixture. In this case, measurement of the mean 
square size in any one direction determines the total energy. From Eq. [6j we see 
that by simply measuring the axial mean square size (z 2 )s40 at 840 G and measur- 
ing the axial trap frequency by parametric resonance, we actually measure Z?840, 
the total energy per particle of the strongly interacting Fermi gas at 840 G. This 
determines the total energy per particle in a model-independent way—. 



2.2 Entropy Measurement 

The entropy S of the strongly interacting gas at 840 G is determined by adiabat- 
ically sweeping the bias magnetic field from 840 G to 1200 G, where the gas is 
weakly interacting 38 . The entropy Sw of the weakly interacting gas is essentially 
the entropy of an ideal Fermi gas in a harmonic trap, which can be calculated in 
terms of the mean square axial cloud size (z 2 )i200 measured after the sweep. Since 
the sweep is adiabatic, we have 




(5) 




(6) 



S = Sw ■ 



(7) 



The adiabaticity of the magnetic field sweep is verified by employing a round- 
trip-sweep: The mean square size of the cloud at 840 G after a round-trip-sweep 
lasting 2s is found to be within 3% of mean square size of a cloud that remains at 
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840 G for a hold time of 2s. The nearly unchanged atom number and mean square 
size proves the sweep does not cause any significant atom loss or heating, which 
ensures entropy conservation for the sweep. The background heating rate is the 
same with and without the sweep and increases the mean square size by about 
2% over 2s. The mean square size data are corrected by subtracting the increase 
arising from background heating over the 1 s sweep time^. 

At 1200 G in our shallow trap, we have kpa$ = —0.75, where the Fermi 
wavevector kp = (2mEp/h 2 ) 1 / 2 and the s-wave scattering length as = —2900 
bohr 41 . We find that the gas is weakly interacting: For the lowest temperatures 
attained in our experiments, the gas at 1200 G is a normal fluid that we ob- 
serve to expand ballistically. We have calculated the ground state mean square 
size at 1200 G in our gaussian trap, based on a mean-field theory, (z 2 )wo = 
0.69z^(1200G) 38 , which is close to that of an ideal harmonically trapped gas, 
(z 2 )/0 = 0.75z^(1200G). Here, Zp{B) is the mean square size corresponding to 
the Fermi energy of an ideal noninteracting Fermi gas at magnetic field B, which 
includes the magnetic field dependence of the axial trapping frequency: Ep(B) = 
3m(0 2 (B)z 2 F {B). 

We expect that the entropy of the gas at 1200 G is close to that of an ideal gas, 
except for a mean field shift of the energy. We therefore assume that a reasonable 
approximation to the entropy is that of an ideal Fermi gas, Si((z 2 )i — (z 2 )io), where 
(z 2 )io ls the ground state mean square size of an ideal Fermi gas in the gaussian 
trapping potential of Eq. [3j Here, we apply an elementary calculation based on 
integrating the density of states for the gaussian trap with the entropy per orbital 
^ = -ke[f In f + (1 - /) ln(l - /)], where /(e) is the ideal Fermi gas occupation 
number at temperature T for an orbital of energy e. By calculating Si as a function 
of the difference between the finite temperature and ground state mean square 
cloud sizes, we reduce the error arising from the mean field shift at 1200 G, and 
ensure that 5/ = for the ground state. 

The exact entropy of a weakly interacting gas Sw at 1200G, Sw{(z 2 )w — 
(z 2 )wo), has been calculated using many-body theories 4 ' 7 for the gaussian poten- 
tial of Eq.[3] In the experiments, we determine the value of (z 2 )w — (z 2 )wq, where 
we take (z 2 )wo = 0.7 \zp{ 1200), the value measured at our lowest energy at 1200 
G by extrapolation to T = using the Sommerfeld expansion for the spatial profile 
of an ideal gas. This result is close to the theoretical value, 0.69z^-(1200). 

The entropy versus cloud size curve for an ideal noninteracting Fermi gas and 
the exact value for a weakly interacting gas Sw at 1200 G are plotted in Fig.Q] We 
find that the entropies Sw(A(z 2 ) = (z 2 )w - (z 2 )wo) and Si(A(z 2 ))), agree within 
a few percent over most of the energy range we studied, except at the point of 
lowest measured energy, where they differ by 10%. The results show clearly that 
the shape of the entropy curve of a weakly interacting Fermi gas is nearly identical 
to that of an ideal gas when the mean field shift of the ground state size is included 
by referring the mean square cloud size to that of the ground state. So we have to 
a good approximation, 

S = S w ((z 2 )w - (z 2 )wo) ^ 5/((z 2 ) w - (z 2 )wo). (8) 

Since the corrections to ideal gas behavior are small, the determination of S 1200 
by measuring the axial mean square size (z 2 )i200 relative to the ground state pro- 
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Fig. 1 Color online. Comparison of the entropy versus mean square size curves of a weakly 
interacting Fermi gas with kpas = —0.75 in a gaussian trap with that of a noninteracting Fermi 
gas. Many-body predictions are shown as the green dashed curv&i and the blue dot-dashed 
curve 7 . The red solid curve is our calculation for a noninteracting gas 47 . All of the calculations 
employ the gaussian potential of Eq.[3]with a trap depth Uo/Ep = 10. The points where S = 
automatically coincide, since the entropy is determined in all cases as a function of mean square 
size relative to that of the ground state. This corrects approximately for the mean-field shift of 
the mean square size. Note that for a weakly interacting gas, {z 2 )o is the calculated ground cloud 
size for each theory that includes the mean field energy shift, while for the noninteracting case 
(z 2 )o is (z 2 ) 10, tne value for an ideal trapped Fermi gas. z F is the mean square size for an energy 
E = E F . 

vides an essentially model-independent estimate of the entropy of the strongly 
interacting gas. 



2.3 Sound Velocity Measurement 

Sound velocity measurements have been implemented for Fermi gases that are 
nearly in the ground state, from the molecular BEC regime to the weakly inter- 
acting Fermi gas regime 22 . A sound wave is excited in the sample by using a thin 
slice of green light that bisects the cigar-shaped cloud. The green light at 532 nm 
is blue detuned from the 671 nm transition in lithium, creating a knife that locally 
repels the atoms. The laser knife is pulsed on for 280 /is, much shorter than typi- 
cal sound propagation times ~ 10 ms and excites a ripple in the density consisting 
of low density valleys and high density peaks. After excitation, the density ripple 
propagates outward along the axial direction z. After a variable amount of prop- 
agation time, we release the cloud, let it expand, and image destructively. In the 
strongly interacting regime, we use zero-temperature Thomas-Fermi profiles for 
a non-interacting Fermi gas to fit the density profiles, and locate the positions of 
the density valley and peak. By recording the position of the density ripple versus 
the propagation time, the sound velocity is determined. A detailed discussion of 
potential sources of systematic error is given by Joseph et al.— . 
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For a strongly interacting Fermi gas in the unitary limit, the sound velocity Co 
at the trap center for the ground state is determined by the Fermi velocity of an 
ideal gas at the trap center, Vp = (lEp/m) 1 / 2 and the universal constant j3, 

vf V5 

A precision measurement of the sound speed therefore enables a determination of 
j3— . As discussed below, the values of /3 determined from the E (S) data and the 
sound velocity data are in very good agreement. 



3 Comparison of Thermodynamic Data with Theory 

In the experiments, the raw data consists of the measured mean square cloud sizes 
at 840 G and after an adiabatic sweep of the magnetic field to 1200 G. Using this 
data, we determine both the energy and the entropy of the strongly interacting gas. 
The data is then compared to several recent predictions. 



3.1 Mean Square Cloud Sizes 

We begin by determining the axial mean square cloud sizes at 840 G and after 
the adiabatic sweep to 1200 G. Since the atom number can vary between different 
runs by up to 20%, it is important to make the comparison independent of the atom 
number and trap parameters. For this purpose, the mean square sizes are given in 
units of Zf{B), as defined above. 

The measured mean square sizes are listed in Table Q] In the experiments, 
evaporative cooling is used to produce an atom cloud near the ground state. En- 
ergy is controllably added by releasing the cloud and then recapturing it after a 
short time t] reat as described previously 38 . For a series of different values of th eat , 
the energy at 840 G is directly measured from the axial cloud size according to 
Eq. [6] Then the same sequence is repeated, but the cloud size is measured after an 
adiabatic sweep to 1200 G. In each case, the systematic increase in mean square 
size arising from background heating rate is determined and subtracted. The total 
data comprise about 900 individual measurements of the cloud size at 840 G and 
900 similar measurements of the cloud size after a sweep to 1200 G. To estimate 
the measurement error, we split the energy scale at 840 G into bins with a width of 
AE = 0.04Ep. Measured data points within the width of the energy bin are used 
to calculate the average measured values of the cloud sizes and the corresponding 
standard deviation at both 840 G and 1200 G. 

The ratio of the mean square axial cloud size at 1200 G (measured after the 
sweep) to that at 840 G (measured prior to the sweep) is plotted in Fig.|2]as a func- 
tion of the energy of a strongly interacting gas at 840 G. The ratio is > 1, since for 
an adiabatic sweep of the magnetic field from the strongly interacting regime to 
the weakly interacting regime, the total entropy in the system is conserved but the 
energy increases: The strongly interacting gas is more attractive than the weakly 
interacting gas. A similar method was used to measure the potential energy change 
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<z 2 )84o/4(840) 


(z 2 ) 1200/^(1200) 


Emo/Ef 


5"l200/^iJ 




^1 200 As 


1 


0.568(4) 


0.743(6) 


0.548(4) 


0.63(8) 


0.91(23) 


0.97(5) 


2 


0.612(5) 


0.776(13) 


0.589(5) 


0.99(11) 


1.18(22) 


1.24(9) 


3 


0.661(5) 


0.803(11) 


0.634(5) 


1.22(8) 


1.36(20) 


1.42(7) 


4 


0.697(9) 


0.814(15) 


0.667(8) 


1.30(10) 


1.43(18) 


1.49(8) 


5 


0.74(1) 


0.87(4) 


0.71(1) 


1.6(2) 


1.72(18) 


1.8(2) 


6 


0.79(1) 


0.89(2) 


0.75(1) 


1.7(1) 


1.79(15) 


1-9(1) 


7 


0.83(1) 


0.94(2) 


0.79(1) 


2.0(1) 


2.03(16) 


2.1(1) 


8 


0.89(2) 


1.02(2) 


0.84(2) 


2.3(1) 


2.32(18) 


2.4(1) 


9 


0.91(1) 


1.02(3) 


0.86(1) 


2.3(1) 


2.31(16) 


2.4(1) 


10 


0.97(1) 


1.10(1) 


0.91(1) 


2.55(4) 


2.57(17) 


2.64(4) 


11 


1.01(1) 


1.17(2) 


0.94(1) 


2.74(7) 


2.75(19) 


2.82(7) 


12 


1.05(1) 


1.18(1) 


0.98(1) 


2.78(4) 


2.80(17) 


2.87(4) 


13 


1.10(1) 


1.22(1) 


1.03(1) 


2.89(2) 


2.90(15) 


2.97(2) 


14 


1.25(2) 


1.35(5) 


1.15(1) 


3.21(12) 


3.20(14) 


3.27(12) 


15 


1.28(1) 


1.39(3) 


1.18(1) 


3.28(6) 


3.28(14) 


3.35(6) 


16 


1.44(2) 


1.49(2) 


1.31(2) 


3.49(4) 


3.48(9) 


3.55(4) 


17 


1.53(2) 


1.62(6) 


1.39(2) 


3.74(10) 


3.73(11) 


3.80(10) 


18 


1.58(1) 


1.63(1) 


1.42(1) 


3.76(2) 


3.74(8) 


3.81(2) 


19 


1.70(2) 


1.73(6) 


1.52(1) 


3.94(9) 


3.92(7) 


3.99(9) 


20 


1.83(5) 


1.79(2) 


1.62(4) 


4.03(3) 


4.01(1) 


4.08(3) 


21 


1.93(3) 


1.96(3) 


1.70(2) 


4.28(4) 


4.26(6) 


4.32(4) 


22 


2.11(5) 


2.17(3) 


1.83(4) 


4.55(3) 


4.53(7) 


4.59(3) 



Table 1 Mean square axial cloud size, energy, and entropy measured in a trapped strongly 
interacting Fermi gas. (z 2 )s4o/z 2 r(840) is the scaled axial mean square size at 840 G. 
(z 2 )i20o/Zf(1200) is the scaled axial mean square size at 1200 G. E^o/Ef is the energy per 
particle of a strongly interacting Fermi gas at 840 G, calculated using Eq.|6] Snoo/kB is the cor- 
responding entropy per particle of the gas after an adiabatic sweep of magnetic field from 840 G 
to 1200 G, where the noninteracting entropy curve (for the gaussian trap) is used to determine the 
entropy at 1200 G and the ground state mean square size is assumed to be (z 2 )o = 0.71 4(1200). 
^1200/^fi I s me ideal gas entropy result assuming {z 2 )o = 0.69z 2 r (1200). S* l2 oo/kB is the entropy 
obtained using an exact many-body calculation for kpa = —0.75—. 
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Fig. 2 Color online. The ratio of the mean square cloud size at 1200 G, {z 2 } 1200, to that at 840 G, 
{z 2 )%40, for an isentropic magnetic field sweep. E is the total energy per particle of the strongly 
interacting gas at 840 G and Ef is the ideal gas Fermi energy at 840 G. The ratio converges to 
unity at high energy as expected (the dashed line). Here, zwf is evaluated at 1200 G, while zf is 
evaluated at 840 G (see text). 



in a Fermi gas of 40 K, where the bias magnetic field was adiabatically swept be- 
tween the strongly interacting regime at the Feshbach resonance and a noninter- 
acting regime above resonance 48 . The resulting potential energy ratios are given 
as a function of the temperature of the noninteracting gas^. In contrast, by ex- 
ploiting the virial theorem which holds for the unitary gas, we determine both the 
energy and entropy of the strongly interacting gas, as described below. 



3.2 Energy versus Entropy 



Fig. |3] shows the entropy which is obtained from the mean square size at 1200 G 
(z 2 )i20oAf(1200) as listed in TableQ] First, we find the mean square size relative 
to that of the ground state, ((z 2 )i200- (z 2 )o)Af(1200). We use the measured value 
(z 2 )o = 0.71z|(1200) for the lowest energy state that we obtained at 1200 G, as 
determined by extrapolation to T = using a Sommerfeld expansion for the spatial 
profile of an ideal Fermi gas. Then we determine the entropy in the noninteracting 
ideal Fermi gas approximation: Si[(z 2 )noo~ (z 2 )o\, where we have replaced (z 2 )io 
by the ground state value at 1200 G. As discussed above, this method automati- 
cally ensures that 5 = corresponds to the measured ground state (z 2 )o at 1200 G, 
and compensates for the mean field shift between the measured (z 2 )o for a weakly 
interacting Fermi gas and that calculated (z 2 )io = 0.77 z 2 F for an ideal Fermi gas 
in our gaussian trapping potential. As shown in Fig. Q] the entropy obtained from 
a more precise many-body calculations is in close agreement with the ideal gas 
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Fig. 3 Color online. The conversion of the mean square size at 1200 G to the entropy. The 
dashed line is the calculated entropy for a noninteracting Fermi gas in the Gaussian trap with 
Uo/Ef = 10. (z 2 )o = 0.71 is the measured ground state size for a weakly interacting Fermi 
gas. The calculated error bars of the entropy are determined from the measured error bars of the 
cloud size at 1200 G. Here, zwf is evaluated at 1200 G. 



entropy calculated in the ideal gas approximation. The energy is determined from 
the cloud profiles at 840 G using Eq.|6] 

Finally, we generate the energy-entropy curve for a strongly interacting Fermi 
gas, as shown in Fig. [4] Here, the energy E measured from the mean square axial 
cloud size at 840 G is plotted as a function of the entropy S measured at 1200 
G after an adiabatic sweep of the magnetic field. We note that above S = 4kg 
(E = l.5Ep) the E(S) data (blue dots) for the strongly interacting gas appear to 
merge smoothly to the ideal gas curve (dashed green). 

In addition to the entropy calculated in the ideal gas approximation, Table Q] 
also provides a more precise entropy S* (1200) versus the axial mean square cloud 
size. These results are obtained by Hu et al.— using a many-body calculation for 
kpa = —0.75 at 1200 G in the gaussian trap of Eq.0 



3.3 Testing Predictions from Many-Body Theories 

Perhaps the most important application of the energy-entropy measurements is to 
test strong coupling many-body theories and simulations. Since the energy and 
entropy are obtained in absolute units without invoking any specific theoretical 
model, the data can be used to distinguish recent predictions for a trapped strongly 
interacting Fermi gas. 

Fig. [5] shows how four different predictions compare to the measured energy 
and entropy data. These include a pseudogap theory 8,49 , a combined Luttinger- 
Ward-De Dominicis-Martin (LW-DDM) variational formalism 6 , a T-matrix cal- 
culation using a modified Nozieres and Schmitt-Rink (NSR) approximation 4 -^, 
and a quantum Monte Carlo simulations^. The most significant deviations ap- 
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Fig. 4 Color online. Measured total energy per particle of a strongly interacting Fermi gas at 840 
G versus its entropy per particle. For comparison, the dot-dashed green curve shows E(S) for an 
ideal Fermi gas. For this figure, the ideal gas approximation to the entropy is used, Si20oAfi of 
Tabled 



pear to occur near the ground state, where the precise determination of the energy 
seems most difficult. The pseudogap theory predicts a ground state energy that is 
above the measured value while the prediction of Ref . 6 is somewhat low com- 
pared to the measurement. All of the different theories appear to converge at the 
higher energies. 



4 Determining the Temperature and Critical Parameters 

The temperature T is determined from the measured E(S) data using the funda- 
mental relation, T = dE/dS. To implement this method, we need to parameterize 
the data to obtain a smooth differentiable curve. 

At low temperatures, one expects the energy to increase from the ground 
state according to a power law in T and a corresponding power law in S, i.e., 
E = Eo+aS b . For a harmonically trapped ideal Fermi gas, we have in the Som- 
merfeld approximation an energy per particle in units of Ep given by E = 3/4 + 
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Fig. 5 Color online. Comparison of the experimental energy versus entropy data with the calcu- 
lations from strong coupling many-body theories. Top: The green dashed curve is a pseudogap 
theory^. The black solid curve is a LW-DDM variational calculation 6 . For this figure, the en- 
tropy is given by S\ 2Q0 /kB in Table [T] Bottom: The blue dashed curve is an NSR calculation^ 5 .. 
The red solid curve is a quantum Monte Carlo simulation 7,50 . For this figure, the entropy is given 
by S* noo /k B in TableQ] 



K 2 ( T/Tf ) 2 /2. The corresponding entropy per particle in units of £g is S = 7T 2 ( 7/ 7> ) , 
so that E = 3/4 + S 2 /(2n 2 ) ~ 0.75 + 0.05 S 2 . 



4. 1 Power Law Fit and Temperature for an Ideal Fermi Gas 

We attempt to use a single power law to fit the Ei(S) curve for a noninteract- 
ing Fermi gas in a gaussian trapping potential, with Uo/Ep — 10, as in our ex- 
periments. The energy and entropy are calculated in the energy range 0.75 Ep < 



15 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 



T,/T F 

Fig. 6 Color online. The inset top left shows a single power-law fit to the calculated en- 
ergy versus entropy per particle of a noninteracting Fermi gas in a gaussian trap. The ex- 
act result (blue dots) is well fit by the single power law function (green solid curve) Ej = 
(0.747 ±0.001) + (0.0419 ± 0.0004) S( 2197±0006 '. The inset bottom right shows a two-power 
law fit to the calculated energy versus entropy per particle of a noninteracting Fermi gas in a 
gaussian trap. The exact E/(S) (blue dots), is very well fit by the two power law function (red 
solid curve) of Eq.[l2]for e ^ with a = 0.0480(2), b = 2.08(3), S c = 2.16(3), d = 2.483(7) 
and Eq = 0.7415(2). The temperature of a noninteracting Fermi gas obtained from the fits to 
the calculated energy versus entropy is compared to the exact temperature. The blue dotted line 
corresponds to exact agreement. The green dot-dashed curve is from the single power law fit 
showing deviation only at the highest and lowest temperatures. The red curve is from two power 
law fit of Eq. [12] with e =/= showing excellent agreement over the full range. 



E < 2Ep and displayed as dots in Fig. [6j We find that a single power law Ej — 
(0.747±0.001) + (0.0419 ±0.0004) s(2-i97±0.006) fits the curye yery well oyer this 
energy range. Note that the power law exponent is b ~ 2, close to the low temper- 
ature value. 

Using the fit function, we can extract the reduced temperature Tifn/Tp — 
dEi/dSj as a function of Si and compare it to the theoretical reduced tempera- 
ture T /Tf at the same Si. The results are shown as the green dashed line in Fig. [6] 
We see that the agreement is quite good except below T /Tp =0.1 and above 0.5, 
where the deviation is ~ 10%. 

To improve the fit and to make a more precise determination of the tempera- 
ture, we employ a fit function comprising two power laws that are joined at a cer- 
tain entropy S c , which gives the best fit. When used to fit the data for the strongly 
interacting Fermi gas, we consider two types of fits that incorporate either a jump 
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in heat capacity or a continuous heat capacity at S c . In this way, we are able deter- 
mine the sensitivity of the temperature and critical parameters to the form of the 
fit function. The two types of fits yield nearly identical temperatures, but different 
values of S c and hence of the critical parameters, as discussed below. 

We take the energy per particle E in units of Ep to be given in terms of the 
entropy per particle in units of kg in the form 



£<(S) =E Q + aS b ; 0<S<S c 

E>(S) = Ei + cS d +e(S-S c ) 2 ; S>S C . 



(10) 



We constrain the values of E\ and c by demanding that energy and temperature 
be continuous at the joining point S c : 



£<(5 C ) =E > {S C ) 
dE< \ ( dE > 



dS 



S c 



dS 



(11) 



By construction, the value of e does not affect these constraints and is chosen in 
one of two ways. Fixing e = 0, the fit incorporates a heat capacity jump at S c , 
which arises from the change in the power law exponents at S c . Alternatively, 
we choose e so that the second derivative E"(S) is continuous at S c , making 
the heat capacity continuous. The final fit function has 5 independent parameters 
Eo,S c ,a,b,d, and takes the form 



E<(S) = E +aS b ; 0<S<S, 
E>(S) =E + aS l c 



, b b ( s 



-e(S-S c ) 2 ;S>S c . 



Here, when e is not constrained to be zero, it is given by 

ab 

T 



e= a 4(b-d)S'r 2 



(12) 



(13) 



Fig.[6]shows the improved fit to the calculated energy versus entropy of a noninter- 
acting Fermi gas in a gaussian trap for Uq/Ef = 10, using Eq. 1 121 with e 0, since 
the ideal gas has no heat capacity jump. In this case, both power law exponents b 
and d are close to 2 as for the single power law fit. The temperature determined 
from the fit agrees very closely with the exact temperature, as shown in Fig.|6](red 
solid line). 



4.2 Power Law Fit and Temperature of a Strongly Interacting Fermi Gas 

In contrast to the noninteracting case, we have found that the energy-entropy data 
of a strongly interacting Fermi gas is not well fit by a single power law function—. 
However, the two power-law function fits quite well, with a factor of two smaller 
value of x 2 than for the single power-law fit. Here, we use % 2 = L^^r 1 ) 2 - where 

y (yi) is the fitted (data) value for the i point, and a, is the corresponding the 
standard error. 
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Fig. 7 Color online. To determine the temperature, the energy-entropy data are parameterized by 
joining two power-law fit functions. The red solid line shows the fit that includes a heat capacity 
jump, while the blue dashed curve shows the fit for a continuous heat capacity (see § |4j. For 
comparison, the dot-dashed green curve shows E(S) for an ideal Fermi gas. For this figure, the 
ideal gas approximation to the entropy is used, S1200AB of Table[T] 

Motivated by the good fits of the two power-law function to the ideal gas en- 
ergy versus entropy curve and the good agreement between the fitted and exact 
temperature, we apply the two power-law fit function to the data for the strongly 
interacting Fermi gas. 

Fig. |7] shows the fit (red solid curve) obtained with a heat capacity jump us- 
ing Eq.QHwith e = and a = 0.12(1), b= 1.35(11), d = 2.76(12), the ground 
state energy Eq = 0.48(1), and the critical entropy S c = 2.2(1). Also shown is the 
fit (blue dashed curve) with continuous heat capacity (e 7^ 0) and a = 0.12(2), 
b = 1.31(17), d = 2.9(2), the ground state energy E = 0.48(2), and the critical 
entropy S c = 1.57(29). 

4.3 Estimating the Critical Parameters 

The fit functions for the E(S) data for the strongly interacting Fermi gas exhibit 
a significant change in the scaling of E with S below and above S c . The dramatic 
change in the power law exponents for the strongly interacting gas suggests a 
transition in the thermodynamic properties. The power law exponent is 2.9 above 
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S c (k B ) 


Ec(E F ) 


Tc(T F ) 


Expt£(S) Fit 1 


2.2(1) 


0.83(2) 


0.21(1) 


Expt£(S) Fit 2 


1.6(3) 


0.70(5) 


0.185(15) 


Heat Capacity Experiment* 




0.85 


0.20 


Theory Ref— 






0.30 


Theory Ref Ji 






0.31 


Theory Ref .i 






0.27 


Theory Ref— 






0.29 


Theory Ref— 


2.15 


0.82 


0.27 


Theory Ref^ 


1.61(5) 


0.667(10) 


0.214(7) 



Table 2 Critical parameters for a strongly interacting Fermi gas. The experimental results are 
obtained from fits to the energy versus entropy data using Eq.[l2] Fit 1 uses e = 0, and has a jump 
in heat capacity. Fit 2 constrains e so that there is no jump in heat capacity. The theoretical re- 
sults are presented for comparison. * U sing the present experimental calibration of the measured 
empirical transition temperature, see § 15.31 



S c , comparable to that obtained for the ideal gas, where d = 2.5. The power law 
exponent below S c is 1.35, which corresponds to the low temperature dependence 
E —Eq<* 7' 3 - 8 6 > c i ose t ma t obtained in measurements of the heat capacity, where 
the observed power law was 3.73 after the model-dependent calibration of the 
empirical temperature 37 , see § 15.31 

If we interpret S c as the critical entropy for a superfiuid-normal fluid transition 
in the strongly interacting Fermi gas, then we can estimate the critical energy E c 
and the critical temperature T c = (dE < (S)/dS)s c - For the fits of Eq.[l2]with a heat 
capacity jump (e = 0) or with continuous heat capacity (e 7^ 0), we obtain 

E c =E + aS h c 

T c =abS b c ~ l . (14) 

Using the fit parameters in Eq. [14] yields critical parameters of the strongly inter- 
acting Fermi gas, which are summarized in Table [2] The statistical error estimates 
are from the fit, and do not include systematic errors arising from the form of the 
fit function. 

We note that the fit function for S(E) previously used in Ref. 38 to determine 
the temperature was continuous in S and E, but intentionally ignored the continu- 
ous temperature constraint in order to determine the entropy as a power of E — Eq 
both above and below the joining energy E c . As the continuous temperature con- 
straint is a physical requirement, we consider the present estimate of the temper- 
ature T to be more useful for temperature calibrations and for characterizing the 
physical properties of the gas than the estimate of Ref.—. 

In contrast to the temperature T, the estimate of T c depends on the value of 
the joining entropy S c that optimizes the fit and is more sensitive to the form of fit 
function than the temperature that is determined from the E and S data. For the fit 
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function S(E) used in Ref. , the temperatures determined by the fit function just 
above E c , T c> , and below E c , T c< , were different. An average of the slopes 1 /T c> 
and 1 /T c< was used to estimate the critical temperature. From those fits, the criti- 
cal energy was found to be E c /Ep = 0.94 ± 0.05, the critical entropy per particle 
was S c = 2.7(±0.2)fcg. The estimated critical temperature obtained from the av- 
erage was T c /Tf = 0.29(+0.03/ — 0.02), significantly higher than than the value 
T c /Tp = 0.21 obtained using Eq.Q~2] which incorporates continuous temperature. 

We are able to substantiate the critical temperature T c /Tp = 0.21 by using our 
data to experimentally calibrate the temperature scales in two other experiments. 
In § 15.21 we find that this value is in very good agreement with the estimate we 
obtain by calibrating the ideal gas temperature observed for the onset of pair con- 
densation. Nearly the same transition temperature is obtained in § 15.31 by using 
the E(S) data to calibrate the empirical transition temperature measured in heat 
capacity experiments^. 

Table |2] compares the critical parameters estimated from the power-law fits to 
the E(S) data with the predictions for a trapped unitary Fermi gas from several 
theoretical groups. We note that calculations for a uniform strongly interacting 
Fermi gas at unitarity 51 yield a lower critical temperature, T c /Tp(n) = 0.152(7), 
than that of the trapped gas, where Tp (n) is the Fermi temperature corresponding 
to the uniform density n. Extrapolation of the uniform gas critical temperature to 
that of the trapped gas shows that the results are consistent—. 

Using the parameters from the fits and Eq.[T2l the temperature of the strongly 
interacting Fermi gas, in units of Tp can be determined as a function of the entropy 
per particle, in units of kg, 

T< ( S) = Tc {j~) ;0 ^ 5 ^ 

T > (S) = T c (^-pJ +2e(S-S c );S>S c . (15) 

Here S c is given in Table [2] from the fits to the E(S) data for the strongly inter- 
acting gas, Eq. [14] gives T c . Fig. [8] shows the temperature as a function of entropy 
according to Eq. Q3] for fits with a heat capacity jump and for continuous heat 
capacity. 



5 Temperature Calibrations 

The estimates of the temperature of the strongly interacting Fermi gas as a function 
of the entropy can be used to experimentally calibrate the temperatures measured 
in other experiments, without invoking any specific theoretical models. The JILA 
group measures the pair condensate fraction in a strongly interacting Fermi gas of 
40 K as a function of the initial temperature 7} c in the noninteracting regime above 
the Feshbach resonance 26,31 . In these experiments, a downward adiabatic sweep of 
the bias magnetic field to resonance produces a strongly interacting sample. Using 
our E(S) data, we relate the endpoint temperatures for adiabatic sweeps of the bias 
magnetic field between the ideal and strongly interacting Fermi gas regimes. We 
therefore obtain the critical temperature for the onset of pair condensation in the 
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Fig. 8 Color online. The temperature versus the entropy of a strongly interacting Fermi gas from 
the fits to the E(S) measurement. The red solid curve is given by Eq.[l5]for e = (heat capacity 
jump) and the blue dotted curve is for e ^ (continuous heat capacity). 

strongly interacting Fermi gas, and find very good agreement with our estimates 
based on entropy-energy measurement. 

In addition, we calibrate the empirical temperature based on the cloud profiles, 
which was employed in our previous measurements of the heat capacity—. 

5.1 Endpoint Temperature Calibration for Adiabatic Sweeps 

We relate the endpoint temperatures for an adiabatic sweep between the strongly 
interacting and ideal Fermi gas regimes. Eq.[l5]gives the temperature of the strongly 
interacting gas as a function of entropy, i.e., T(S). 

Next, we calculate the entropy per particle 5/(7}) for an ideal Fermi gas in our 
gaussian trap, in units of kg , with 7} in units of Tp, as used in § 14. 1 I to determine 
7i/(5/). For an adiabatic sweep between the strongly interacting and ideal Fermi 
gas regimes, where 5 = 5/, the temperature of the strongly interacting gas is related 
to that for the ideal Fermi gas by 



T = T[S I {T I )], 

which is shown in Fig. [9] 



(16) 



21 




Fig. 9 Color online. Experimental temperature calibrations: Temperature of a strongly interact- 
ing Fermi gas (T/Tp) compared to the temperature of a noninteracting Fermi gas (Tj/Tf) for an 
adiabatic sweep between and strongly interacting and ideal Fermi gas regimes (equal entropies). 
The solid red curve is obtained from the fit to the E(S) data with a heat capacity jump and the 
dashed blue curve is obtained from the fit with continuous heat capacity. The dashed green curve 
shows the value of the empirical temperature 7y„ , as obtained from the cloud profiles in Ref— , 
versus the corresponding reduced temperature of the strongly interacting Fermi gas at the same 
energy. The dotted black line denotes equal temperatures. 



For an adiabatic sweep from the ideal Fermi gas regime to the strongly inter- 
acting Fermi gas regime at low temperature T <T C , the reduced temperature of 
the strongly interacting gas is greater than or equal to that of the ideal gas. This 
arises because the entropy of the strongly interacting gas scales as a higher power 
of the temperature than that of the ideal gas. 

In our present experiments, we could not take data at high enough temperatures 
to properly characterize the approach of the temperature to the ideal gas regime. 
Above T c , our E(S) data are obtained over a limited range of energies E < 2Ep to 
avoid evaporation in our shallow trap. In this energy range, our data are reasonably 
well fit by a single power law. However, such a power law fit cannot completely 
describe the higher temperature regime. We expect that the temperatures of the 
strongly interacting gas and ideal gas must start to merge in the region S = 4kg, 
where the E(S) data for the strongly interacting gas nearly overlaps with the E(S) 
curve for an ideal gas, as shown in Fig.lH 
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From Fig. [8] S > 4kg corresponds to T /Tp > 0.6, approximately the place 
where the calibrations from the two different power law fits (for e = and e ^ 0) 
begin to differ in Fig. [9] We therefore expect that the single power law fit overes- 
timates the temperature T of the strongly interacting gas for T > 0.6 Tp, yielding 
a trend away from ideal gas temperature, in contrast to the expected merging at 
high temperature. 



5.2 Critical Temperature for the Onset of Pair Condensation 

In Ref. 26 , projection experiments measure the ideal Fermi gas temperature 7) c 
where pair condensation first appears. In those experiments, 7} c is estimated to 
be 0.18(2), 7> 31 , From the calibration, Fig. [9] we see that for T Ic = 0.187>, the 
corresponding temperature of the strongly interacting gas is T c = 0.19(2)7f for 
both the red solid and blue dashed curves, which is almost the same as the ideal 
gas value. The critical temperature of the strongly interacting gas for the onset 
of pair condensation is then 0.19(2)7^, in very good agreement with the values 
T c = . 2 1 ( 1 ) Tp and T c = . 1 85 ( 1 5 ) 7> that we obtain from the two fits to the E (S) 
measurements. This substantiates the conjecture that the change in the power law 
behavior observed at T c in our experiments corresponds to the superfluid transi- 
tion. 



5.3 Calibrating the Empirical Temperature obtained from the Cloud Profiles 

In our previous study of the heat capacity, we determined an empirical tempera- 
ture Tfi, = T as a function of the total energy of the ga s 37 i 19 . The gas was initially 
cooled close to the ground state and a known energy was added by a release and re- 
capture method. Then a Thomas-Fermi profile for an ideal Fermi gas was fit to the 
low temperature cloud profiles to determine the Fermi radius. Holding the Fermi 
radius constant, the best fit to the cloud profiles at higher temperatures determined 
the effective reduced temperature, which is denoted f. The E (f) data 37 was ob- 
served to scale as E — Eq = 1.54Ef T 143 for T > 0.33, while below T = 0.33, 
the energy was found to scale as E — Eq — 4.9 Ep f 2 53 . The transition point oc- 
curs at an energy E c = 0.85 Ep, which is close to the value 0.83Ep obtained from 
power-law fit to the E(S) data for the fit with a heat capacity jump. Assuming that 
T c = 0.33 corresponds to the superfluid-normal fluid transition, we can determine 
the corresponding value of T c /Tp for the strongly interacting gas. 

To calibrate the empirical temperature we start with E (f). Then, as discussed 
in § 17.11 Eq. 1 121 determines E(T) and hence T{E)ivova the fits to our E(S) data. 
Hence T ( T ) = T [E ( T )] , where T = T / Tp is the reduced temperature of the strongly 
interacting gas and E = E /Ep is the reduced energy. For simplicity, we give the 
analytic results obtained using the e = fit to the E{S) data, 

— =0.42f a66 ; 0<f <0.33 
T F - - 

^- = (0.80f 143 - 0.09) 064 ; f > 0.33. (17) 
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Fig. [9] shows the full calibration (green dashed curve). For comparison, the cal- 
ibration obtained from the pseudogap theory of the cloud profiles gave T /Tp — 
54 f 0.67 for f < o.33, and T/T F = 0.71 T above f = 0.33. For % = 0.33, we 
obtain from Eq. \\l\ TrlTw = 0.20 (see Fig. |9J), in good agreement with the value 
obtained for the onset of pair condensation and with the values T c = 0.185(15) 
and T c = 0.21 (1) determined from the fits to the E(S) data. 

6 Measuring the Ground State Energy 

Measurement of the ground state energy of a unitary Fermi gas provides a strin- 
gent test of competing many-body theoretical predictions and is therefore of great 
interest. For a unitary Fermi gas of uniform density in a 50-50 mixture of two spin 
states, the ground state energy per particle can be written as 

E g = (l + p)^e F (n), (18) 

where Ep (n) is the local Fermi energy corresponding to the density n. The ground 
state energy of the unitary Fermi gas differs by a universal factor £, = 1 + j3 from 
that of an ideal Fermi gas at the same density. The precise value of S, has been of 
particular interest in the context of neutron matte r 52 ' 53 ' 54 ' 55 , and can be measured 
in unitary Fermi gas experiments^^. 

The sound speed at temperatures near the ground state determines j3 according 
to Eq. [9) We have made precision measurements of the sound speed in a trapped 
Fermi gas at the Feshbach resonance^ 2 -. At 834 G, we vary the density by a factor 
of 30 to demonstrate universal scaling and obtain the value co/vp = 0.362(6). 
Using Eq.|9]then yields /3 = —0.565(15). Note that the reference Fermi velocity 
Vp depends on the Fermi energy of an ideal gas at the trap center and hence on both 
the trap frequencies and atom number (as N 1 ^ 6 ), which are carefully measured to 
minimize systematic errors 22 . While the energy of the gas as measured from the 
mean square cloud size was close to the ground state value, the precise temperature 
of the gas was not determined. 

The universal parameter j5 also can be determined by measuring the ground 
state energy Eq of a harmonically trapped unitary Fermi gas, which is given by 

E Q = ^y/T+PE F . (19) 

Our E(S) data enables a new determination of Eq by extrapolating the measured 
energy E(S) to S = 0. As pointed out by Hui et al.— , this method avoids a system- 
atic error arising when the finite temperature is not determined in the measure- 
ments. 

From both of our fit functions below S c , we obtain Eo/Ep = 0.48(1). Eq. [T9l 
yields f} = —0.59(2). This result is slightly more negative than that obtained in 
the sound speed experiments, which is reasonable since the sound speed measure- 
ments are done at finite temperature. Both results are in very good agreement. 

One possible systematic error in these measurements arises from the determi- 
nation of the atom number. The measurements of /3 from the sound speed and from 
the energy-entropy measurements were done in different laboratories. The close 
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agreement is gratifying, considering that the imaging systems that determine the 
atom number employed a -polarized light for the sound speed experiments, while 
the entropy-energy measurements used x-polarized light, for which the resonant 
optical cross section is a factor of two smaller than for a polarization. To exam- 
ine the systematic error arising from the atom number determination, we employ 
a third method to measure /3 based on the measured ratio of the cloud size at 840 
G and at 1200 G, which is number independent. 

The ratio of the ground state mean square sizes for the weakly and strongly 
interacting gases is predicted to be 



Note that we obtain (z 2 )o,120oAfi200 = 0-^9 from a mean field calculation—, in 
agreement with that obtained using a many-body calculation^. 

Our measurements for the ground state mean square size at 1200 G are ac- 
complished by fitting a Sommerfeld expansion of the axial density for an ideal 
Fermi gas to the cloud profile 38 ' 56 . The fit determines the Fermi radius a z and re- 
duced temperature T/T F , yielding (z 2 ) ,i200 = cr 2 /8 = 0.7l4(1200) for T = 0, 
close to the predicted value of 0.69. The ground state energy E$ = 0.48 Ep (840) 
at 840 G from the entropy-energy experiments determines the ground state mean 
square size as (z 2 )o,840 = 0.48 z^^O). Hence, r Q = 0.71/0.48 = 1.48. The cor- 
responding fi = —0.61(2) from Eq. [20j Since the mean square sizes are deter- 
mined from the images and the ratio z 2 7 (840)/z 2 7(1200) is number independent, 
this result shows that the systematic error arising from the number measurement 
is within the quoted error estimate. 

We also can determine f} by directly extrapolating to zero entropy the ratio of 
the axial mean square size of the weakly interacting Fermi gas at 1200 G to that 
of strongly interacting gas at 840 G. When this is done, we obtain j5 = —0.58, 
in very good agreement with the estimates based on the sound speed and ground 
state energy. 

Finally, we can estimate the correction to the ground state energy arising from 
the finite scattering length at 840 G, as = — 73616ao- For the trap conditions in 
the E(S) measurements, kpcig = —18, where kp is the wavevector for an ideal 
Fermi gas at the trap center. To estimate the true unitary ground state energy 
at a$ = °°, we first determine the leading order l/(agkp (n)) correction to the 
trapped atom density, where kp (n) is the local Fermi wavevector corresponding 
to the density n. The local chemical potential is estimated from Ref. 57 . Using 
the notation of Eq. [3] and a harmonic approximation, the corrected density yields 
(r 2 /a 2 ) = (3/8)[l - (128/105^(0.64/^^)], where a is the Fermi radius for 
the unitary gas. According to the virial theorem (see Eq.[5]), the mean square size 
and energy of the unitary gas are corrected by the same factor. The unitary ground 
state energy is then 



ro = 



(z 2 )o,i20o/z 2 ,(1200) _ 0-69 



(20) 



(z 2 )o,84o/z^(840) (3/4) V / TT 7 3' 



E (k F a s ) 



(21) 



1- (128/105710(0.64/^5)' 



Former, = -18, we obtain £ (°°) =0.986 E (- 18) and the value of j3 = -0.59(2) 
obtained directly from E(S = 0)=E = 0.48(1) is shifted to j5 = -0.60(2). We 
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P 


E(S) Experiment 


-0.59(2) 


Sound Velocity Experiment 


-0.565(15) 


Cloud Size Ratio Experiment 


-0.61(2) 


Ref.5^59 


-0.58(1) 


Refi2 


-0.56(3) 


Ref4 


-0.599 


Ref. 60 


-0.60(1) 


Ref£ 


-0.646(4) 



Table 3 Universal interaction parameter /3 . 



also obtain the corrected value of r = 0.71/(0.986*0.48) = 1.50 in Eq.|20land 

P = -0.62(2). 

Table[3]compares the values of /3 obtained in our experiments to several recent 
predictions. Note that the table does not include the finite kp a correction for the 
E(S) measurement at 840 G described above. 



7 Universal Thermodynamic Functions 

Using the E(S) data for the strongly interacting Fermi gas and the temperature 
determined from the two power-law fits, we estimate several universal functions. 
First, we determine the dependence of the energy on temperature E(T) and the 
corresponding heat capacity, C(T). Then we find the global chemical potential of 
the trapped gas as a function of the energy ju ? (Zs). 



7.1 Energy versus Temperature 



The energy is readily determined as a function of temperature using Eq.[l2]for the 
case where there is a heat capacity jump and e = 0, 



E < (T)=E C + 
E > {T)=E C + 



S C T C 




b 




ScT c 




d 





; < T < T c 

;T>T C1 



(22) 



where the energy (temperature) is given in units of Ep (Tp) and the critical energy 
E c is 

E c =E + S c T c /b, (23) 
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Fig. 10 Color online. The energy of a strongly interacting Fermi gas versus temperature, from 
the fits to the E(S) data. The red curve shows E(T) as determined from the fit with a heat 
capacity jump (e = 0) in Eq.[T2] The blue dashed curve shows E(T) as determined from the fit 
with continuous heat capacity (e ^ 0). 

with Eq the ground state energy. For the case with e ^= 0, where the heat capacity 
is continuous, we determine the ordered pairs [E(S),T(S)] as a function of S and 
plot E(T). Fig.llOlshows the results using the best fits for both cases. 

Of particular interest is the low temperature power law. For e = 0, we obtain 
b = 1.35 and bj(b — 1) = 3.86. Since b is near 4/3, the energy relative to the 
ground state scales approximately as T 4 . This is consistent with sound modes 
dominating the low energy excitations. However, one would expect instead that 
the free fermions on the edges of the trapped cloud would make an important 
contribution to the low energy excitations 8 . Over an extended range of T < T c , 
the net entropy arising from the Bose and Fermi excitations has been predicted to 
scale as T 2 , yielding an energy scaling 8 as E — Eq « T 3 . In this case, one would 
expect that E -E S 3 / 2 , i.e., b = 3/2 in Eq.El so that b/(b- 1) = 3. Hence, 
the low energy power law exponents for the entropy should be between 4/3 and 
3 /2, which is barely distinguishable for our data. 



7.2 Heat Capacity versus Temperature 

The heat capacity at constant trap depth C = dE/dT is readily obtained from 
Eq. [22] (where there is a heat capacity jump, since we have constrained e = in 
Eq. 112b. For this parameterization, 
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Fig. 11 Color online. Heat capacity versus temperature given by Eq.[24]for a strongly interacting 
Fermi gas. The red dot-dashed curve shows the heat capacity when there is a jump at T c /Tp = 
0.21. The blue solid curve shows the heat capacity when the heat capacity is continuous. For 
comparison, the light-blue dashed curve shows the heat capacity obtained for an ideal Fermi gas 
(using the fit function of Fig. [7). 



S ( T \ ~ 

C>{T) = d~l{T c ) (24) 

where T and T c are given in units of Tp, and S c is given in units of kg- For the 
fit with a continuous heat capacity, we use T(S) to find C(S) = T(S)/(dT/dS), 
and plot the ordered pairs [C(S),T(S)]. The heat capacity curves for both cases 
are shown in Fig. II 11 

According to Eq.|24l a jump in heat capacity occurs at S c : C < (T C ) = S c / (b — 
1) and C>(r c ) = S c /(d — 1) differ when the power law exponents b and d are 
different. This is a consequence of the simple two power-law structure assumed 
for the fit function E(S) given by Eq.[T2]for e = 0, and cannot be taken as proof 
of a true heat capacity jump. At present, the precise nature of the behavior near 
the critical temperature cannot be determined from our data, and it remains an 
open question whether the data exhibits a heat capacity jump or a continuous heat 
capacity. 




7.3 Global Chemical Potential versus Energy 

The global chemical potential /j, g is readily determined from the fits to the E(S) 
data for a strongly interacting Fermi gas, which obeys universal thermodynamics. 
The local energy density generally takes the form e = T s + jin — P, where £ is the 
local internal energy, which includes the kinetic energy and the interaction energy. 
Here, n is the local density, fj. is the local chemical potential, P is the pressure and 
5 is the total entropy per unit volume. 
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The local chemical potential can be written as )i = /J. g — U, where U is the 
trap potential. In the universal regime, where the local pressure depends only on 
the local density and temperature, we have P = 2e/3, as noted by Ho 15 . Hence, 
5e/3 = T s + [jx g — U)n. Integrating both sides over the trap volume and using 
/ d 3 xe=NE-N(U), where E and (£/}) are the total energy and average potential 
energy per particle, respectively, we obtain 

jNE - ^N(U) =NTS + ll g N, (25) 

where S is the entropy per particle. For simplicity, we assume harmonic confine- 
ment and use the virial theorem result, (U) =E/2, from Eq.[5] which holds in the 
universal regime. Then, Eq.[25]yields the global chemical potential of a harmoni- 
cally trapped Fermi gas in the universal regime, 



Pg = -E-TS. (26) 

By using the fit to the measured entropy-energy data to obtain the tempera- 
ture T = dE I dS from Eq. [I2j the global chemical potential of a trapped unitary 
Fermi gas can be calculated from Eq.|26] For e = 0, where the heat capacity has a 
jump, the simple power law fits above and below E c each yield a different linear 
dependence of fi g on E, 

H g (E) = \e + ( \ - b ) (E - Eq); E < E < E c 



pi g {E) =ii g {E c )+[~-d\{E-E c ); E> E c , (27) 



3 V 3 

3 

where n g {E c ) = 4E Q /3 + (4/3 - b) (E e ~ Eq). 

We plot the chemical potential in Fig. [12] The data points are obtained using 
Eq. [26] with the measured energy E and entropy S and the temperature determined 
from the fit to the E(S) data, using e = in Eq.Q~2] i.e., with a heat capacity jump. 
The solid red curve is given by Eq. [27] We note that the low temperature data 
points in E(S) are best fit with the power law b = 1.35, which is close to 4/3. 
According to Eq.[27] this produces a slope near zero for Eq < E < E c . Since the 
power-law fit above E c gives d = 2.76, the slope according to Eq.[27]changes from 
nearly zero for Eo < E < E c to negative for E > E c . 

Note that from Eq.[26] we obtain the slope 

Since the entropy S is continuous, we see that a jump in the heat capacity produces 
a corresponding jump in the slope of /i ? versus E. 

For comparison, Fig. [12] also shows the chemical potential obtained for e 7^ 
in Eq.[T2] where the heat capacity is continuous (blue dashed curve). 



29 




0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 



E/E F 

Fig. 12 Color online. The global chemical potential versus the total energy of a strongly inter- 
acting Fermi gas. The data points are calculated from the measured E — S data and the fitted T, 
where T is determined by the fit parameters in Eq.[l2]for e = 0. The standard deviation for each 
point of the chemical potential is determined by the standard deviation of the measured E — S 
data. The solid red lines (heat capacity jump) and blue dashed curve (heat capacity continuous) 
are determined by the fit parameters used in Eq. [12] according to Eq.[27] The green dot-dashed 
curve shows the ideal Fermi gas result for the fit function of Fig. ??. 



8 Conclusion 



We have studied the thermodynamic properties of a strongly interacting Fermi gas 
by measuring both the energy and the entropy. The model-independent data ob- 
tained in both the superfluid and the normal fluid regimes do not employ any 
specific theoretical calibrations, and therefore can be used as a benchmark to 
test the predictions from many-body theories and simulations. Parameterizing the 
energy-entropy data determines the temperature of the strongly interacting Fermi 
gas and also yields estimates of the critical parameters. We use the measured data 
to calibrate two different temperature scales that were employed in observations 
of the onset of pair condensation and in heat capacity studies. These calibrations 
yield critical temperatures in good agreement with the results estimated from our 
energy-entropy data. Our data does not determine whether the heat capacity ex- 
hibits a jump or is continuous at the critical temperature. However, for a finite 
system with nonuniform density, the latter is most likely. Considering that there 
is huge interest in determining the detailed behavior of the superfluid transition 
in a strongly interacting Fermi gas 61 , more precise determinations of the critical 
temperature, the heat capacity, and the chemical potential near the critical point, 
as well as the high temperature behavior and the approach to the ideal gas limit, 
will be important topics for future research. 
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